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Abstract 

We apply the quark delocalization and color screening model to nucleon- 
baryon scattering. A semi-quantitative fit to N-N, N-A and N-E phase 
shifts and scattering cross sections is obtained without invoking meson ex- 
change. Quarks delocalize reasonably in all of the different flavor chan- 
nels to induce effective nucleon-baryon interactions with both a repulsive 
core and with an intermediate range attraction in the cases expected. 
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I. INTRODUCTION 



QCD has acquired significant experimental support as the correct fundamental theory of 
the strong interaction. However in the low energy region, its nonperturbative nature makes 
it hard to use directly for the study of complicated systems such as hadron interactions 
and exotic quark states. QCD inspired models are generally used in these cases. Meson 
exchange models @D, on the other hand, are quite successful in explaining the hadron 
interactions. But since the quark and gluon degree of freedom are completely integrated 
out and many phenomenological coupling constants and short-range phenomenology are 
involved in this model, it is hard to use it to make predictions of the properties of exotic 
quark states. Hybrid meson-gluon exchange quark cluster models have been developed by 
few groups iHl and have achieved a Quantitative fit to the scattering data that accounts 
well for nucleon-baryon interactions 0. However, in ref. || baryons are assumed to be so 
stiff that no internal distortion would be induced, no matter how close the interacting 
baryons are. The Hamiltonian used is a direct extension of that used in single hadrons. 
It is possible that some physics has been precluded by these assumptions a priori. For 
example, the QCD vacuum in between baryons can be expected to vary as the quark 
matter density increases as two colliding baryons approach each other and quark 
percolation between hadrons may occur at short distances. On the other hand, the color 
confinement interaction is screened at large distances due to excitation of quark-antiquark 
(qq) pairs, as has been shown by unquenched lattice QCD calculation i. Also, various 
kinds of multigluon exchange interactions cannot be included in two body confinement or 
the Fermi-Breit form for the interaction; nor can they be studied in the pure valence quark 
model of single hadrons. An example of such a three gluon exchange interaction, which is 
impossible for a qq meson and does not contribute within a colorless q 3 baryon but does 
contribute to hadron-hadron interactions and to multiquark states, has been discussed in 
^. [0- 

Recently, the proposed d! (IJ P = 00 ) dibaryon has been studied by Faessler's group H, 
who concluded that both bag and hybrid meson-gluon exchange models cannot obtain such 
a state with a mass as low as 2.06GeV. If the d' is proven to exist, that will argue strongly 
against the completeness of these models. Moreover, some well established facts are either 
impossible or hard to address within either the meson exchange or the hybrid meson-gluon 
exchange model approaches. For example, nuclear and molecular forces have been known 
to be similar for more than half a century % . If the nuclear intermediate range attraction is 
due to meson exchange, then the similarity would be accidental because the molecular force 
certainly cannot be due to electron-positron pair exchange. The nucleus as a collection of 
nucleons which has been proven to be a good approximation through the nuclear structure 
studies, which leaves a question as to why nucleus is not simply a collection of quarks. 
A different approach has been developed by our group, which we call the quark 
derealization, color screening model (QDCSM) 113. It is aimed at including more of the 
physics needed for the study of exotic quark states. With the lessons learned in earlier 
studies of such states in mind, we first tested our model on N-N scattering data 0. We 
found that the QDCSM is able to produce both the N-N short range repulsion and 
intermediate range attraction simultaneously without invoking meson exchange. The 
similarity between nuclear and molecular forces thus obtains a natural explanation, i.e., 
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the intermediate range attraction is due to distortion of the internal structure of the 
constituent nucleons (and atoms, respectively). Nuclear collections of Qq, 9q and 12q 
systems can be shown to be energetically favored when they have structures organized into 
something close to the conventional form of d, 3 H, 3 He and 4 He. On the other hand, a 6q 
system with IJ P = 03 + quantum numbers is energetically favored to be in a six quark 
state, d*, instead of two As. This state is predicted to be a narrow resonance with width of 
order MeV ©. 

Only very limited experimental data is available for nucleon-hyperon (N-Y) scattering. 
Nevertheless it provides a further check of hadron interaction models. Here we report the 
results of applying the QDCSM to N-Y scattering. We should emphasize immediately that 
we are not so ambitious as to expect to achieve a quantitative fit to all nucleon-hyperon 
scattering with an almost parameter free model; more modestly, we only seek to test 
further whether the QDCSM contains basically the right physics. 

II. QUARK DELOCALIZATION AND COLOR SCREENING MODEL 

We present a two baryon system as an example to illustrate the QDCSM. The generator 
coordinate method (GCM) is used to describe the Qq system. The GCM basis wave 
function (WF) is assumed to be 



where A is the normalized antisymmetrization operator and Xci) r ls i T i {i — 1,2) are color, 
spin-isospin WFs. [ }st means that SxTx, S2T2 are coupled to channel spin S and isospin T. 
Xa is always the 3q color singlet state. r] S . Ti is the SUq t symmetric spin-isospin WF for 
the N-N channel, but for A and S, we use the uds spin-flavor asymmetric hyperon WF, 



where the arrows refer to the quark and overall baryon spin, and antisymmetrization will 
be applied to the five u, d quarks only for N-A and N-E channels. This choice explicitly 
distinguishes s quarks and allows for flavor symmetry breaking effects to be calculated 
more easily, but the results are the same as those which use an SU% symmetric baryon 
WF and totally antisymmetric six quark states. 
Our spatial WFs take the form 



tf(l---6) 

<M123) 
^ 2 (456) 



^foMl23)V>B a (456)]sr, 

Xci (123)^(123)^(1)^(2)^(3), 

Xci (45% 52 T 2 (456)^(4)^(5)^(6) , 



(1) 
(2) 
(3) 




(4) 



M?) = (M?) + e<Mr)) / V /l + e 2 + 2e(0 L |(/) i? ) 



(5) 
(6) 



Mr) = (M**) + eMr)) I V 1 + e ' + HMr) 
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0,(f)=(— ) (7) 

/ 1 \ 3 / 4 (r-? 2 ) 2 

where b is the size parameter of baryon. s = s 2 — Si, which is the separation of two 
reference centers, plays the role of the generator coordinate in our model. The only 
difference from the usual GCM basis WF is found in the single quark orbital WF eq.(§, ||): 
A derealization parameter e(s), which will be determined variationally by the six quark 
dynamics for each separation s = \s\, has been introduced to describe the mutual 
percolation of quarks originally confined in different baryons. This basis WF includes the 
six quark bag-model-like WF (e = 1 case, which is spherical only for s = 0) and the usual 
quark cluster model WF (e = case) as two extremes. Note that intermediate 
configurations corresponding to mutually distorted baryons are allowed within the 
variational Hilbert space. 

The Hamiltonian of the six quark system is taken as 

H(l ■ ■ ■ 6) = EK + £-) + £ (y$ + V?) , (9) 

i=l i i<j=l 



1 ird(r) / 2^ ^ 
— 1 + o°i ■ °i 

Tij VTLiVTLj \ 6 



(10) 



where A(a) is the SU% Gell-Mann (SU% Pauli) operator, r = r» — fj and the other symbols 
have their usual meaning. Momentum dependent and tensor interactions are neglected in 
this calculation. is the Fermi-Breit approximation to single gluon exchange and the 
color screened confining interaction is defined by 

g A A | r if occur in the same baryon orbit, 

13 c 1 3 \ l ~ e K Kr if i,j occur in different baryon orbits, 

Explicitly, eq. (|TTl) means we use the non-screened, color- confinement potential, namely 
linear confinement, to calculate matrix elements (LL\V\LL), (RR\V\RR), (LL\V\RR) and 
(RR\V\LL) , and use the color-screening confinement potential to calculate the other 
matrix elements, such as (LR\V\LR), (LL\V\LR) . Here 
(LL\V\LL) = (M^)Mj)\V tJ \M^Mj))^tc. 

For ( LR | V | LR), the interacting quarks are clearly always in different baryons and so 
screened at longer distances. For (LL\V\LR) , which form applies is not obvious and a 
model choice must be made. One could consider other choices such as an average of the 
two confinement forms. Another ambiguity occurs for s — 0, where we have only one 
baryon orbit. We define the matrix elements at s = to be the limit of the values as s 
tends to zero. Both of these model choices are consistent. We recognize that making these 
additional assumptions means that we no longer have simply a potential model but rather 
that we are implicitly including an approximation to many-body, low energy QCD 
interactions which cannot be included in two body confinement and Fermi-Breit 
interactions. (A similar inclusion occurs in models such as the "flip-flop" model. 

@) In 
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this sense, we are extending an effective matrix element method from bound states to 
scattering states. 

A good feature of this model Hamiltonian is that it reduces to the usual non-screening, 
color-confinement model Hamiltonian for a single hadron and for two hadrons in the 
asymptotic region but the spurious color van der Waals force have been eliminated. The 
historical triumphs of the constituent quark model in explaining hadron spectroscopy are 
retained and the model parameters m (u,d quark mass), a s (quark-gluon coupling 
constant), a c (strength of the confinement potential) and b (baryon size parameter) can be 
determined by the nucleon mass, N — A mass splitting, and the stability condition for 
nucleon size, dM^(b)/db = 0, along with the usual choice, m = Mjv/3. The s quark mass 
(m s ) is most accurately determined by the difference of A and E masses. 
The color screening constant k is directly taken from lattice QCD results O and this is the 
reason for our having chosen a linear confinement and exponential color screening in this 
calculation even though quadratic confinement may be more proper for a nonrelativistic 
model 0. The parameters fixed in this way are: 

m = 313MeV, m s = 521.7MeV, b = 0.625fm, 

a s = 1.71, a c = 39.1MeV fm -1 , k = l.llllfm -1 (12) 

Note that the A and E masses calculated are 1025 MeV and 1103 MeV, 90 MeV lower 
than experimental values. 



III. CALCULATION METHOD 

Due to derealization, our GCM basis WF, eq.(ffl), includes not only the usual q 3 -q 3 
clustering, but also q 6 , q 5 -q and q 4 -q 2 clustering, and therefore can not be factorized into 
internal, relative and center of mass WF in the interaction region. The usual cluster model 
method has to be extended @. Suppose ^ is a solution of our model Hamiltonian eq.@, 

(H-E)V = 0. (13) 

In general, both local and nonlocal interactions are included in H, but the nonlocal 
interaction is nonzero only within a limited interaction region r < a, where a is roughly 
determined by the scale at which the overlap of different orbitals, Eqs.(^) and (§), becomes 
negligible, and r refers here to the separation between two three-quark clusters (see 
Eq.(pTD below). If we separate the whole space into interaction and asymptotic regions, 
then in the interaction region, we can rewrite eq.([T3Das, 

(H + £- E)^ = C^, (14) 

where C is the Bloch operator @, which was introduced by Bloch to make the 



Hamiltonian Hermitian in a finite space. Taking the Hermitian conjugate of eq. ([14|) and 
using W + £t = h + C, we obtain 

(*\H-E\* t )\%=(*\tf-£\%)\Z. (15) 
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where the notation |g is intended to convey that the integration in r is restricted to 

< r < a and \l/ t is a trial WF. Then a variational functional, J(\I/ t )|o, can be defined in 

the interaction region as, 

J{%)\1 = (%\H - E\%)\% - (*\H - E\%)\ a 

= (%\H - E\V t )\™ - (%\H - E\%)\? - (tf |£* - C\%)\ a , (16) 

and this functional does have a variational minimum with respect to variation of the trial 
WF, ty t . This can be seen from the fact that the first line is quadratic in the difference 
between the trial and exact WFs since, due to eq. (p~3|) , 

(tyt\H ~ E\^)\q = (ty\H — E\^)\q = 0. The second line simply involves rewriting the first 
term in terms of the full and exterior ranges, and the second term has had the substitution 
made from eq . (|15|) . 
The GCM WF is written as 

*? CM = I f(m^)ds, (17) 



where &(s) is the six quark WF of eq.([|). Upon substituting this trial WF in eq. flTB]) and 
doing a partial wave decomposition (Only central interactions are studied here; if 
non-central interactions are to be included, this partial wave decomposition would have to 
be extended correspondingly), we obtain 

J,(**)|g = / dsds'f^f^K^^s') - Ct(a), (18) 

where 

Kf CM {s, a') = K? CM {s, s') - K[ GCM {s, s'), (19) 

K ? CM ^ S,) = f Yr m (s){^(s)\H - EMtyY^ffldstf, (20) 
ss' J 

h 2 ( d 2 1(1 + 1) 



/■oo 

K' l GCM (s,s') = c l drTfas) 

J a 



2/i \dr 2 r 2 



+ V c (r) - E r 



Ti(r,s') (21) 



where here, f = (&±&±a _ n+5+f 3 ^ Cl = \_ (—) s+t+1 5ab, anc j ^ ig ^ e energy of 
relative motion. 

To obtain eq.(|H]), we have assumed that, in the asymptotic region: All of the exchange 
color interactions have died out; only the long range Coulomb interaction, V c (r), for 
charged baryons may remain; the derealization has disappeared (e = 0); and the six quark 
system has clustered into two three quark baryons A and B. 

Ti(r, s) is the Ith partial wave of the relative motion WF, T(f, s), obtained from the WF of 
eq. ([!]) when it is factorized into internal, relative and center of mass parts of two three 
quark clusters (see eq. (|32[) below). 

/ 1 \ 3 / 4 

r < r >M^) (22) 

/ o \ 3/4 o 

rKr ' S) = (2^) e-^^hnrs^U-z^rs). (23) 
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Finally 



h 2 r 

= (ci—) [9'i{a)9 t i{ a ) -9i{a)9?{a) 



9\{r) 



fi(s)Ti(r,s)ds, 



(24) 
(25) 



and gi(r) is the radial part of the relative motion WF of the exact solution \I/ of eg. (fD|) . 
The Bloch operator C = j^j\S(r — a)^r has been used in deriving eq.(|24]). 
Next we follow Canto and Brink E3 , assuming as a boundary condition that the trial and 
exact logarithmic derivatives are equal at the boundary, 



9_M 
9\{o) 



Li 



9iW 



(26) 



This implies that Ci (a) = 0. 

Making use of the stability property of the functional J(\I/*), we vary J{^ 1 ) with respect to 
the trial WF fi(r) to obtain the equation of motion of fi{r). The advantage of the 
Canto-Brink variational method is that one need not solve this equation of motion, but 
that instead one can obtain the logarithmic derivative L\ directly via 



Eu r z( a ; s i)Qi/ r K a > s j 



Q 



2fiK 



Kl 



'J 



K, 



w(r) 



tr c, 

R GCM { 
1(1 + 1] 



dr 



Ti(r, Si)w(r)Ti(r, Sj) 



(27) 
(28) 



■Si, Sj y 



+ ¥ [VC{r) 



E r 



Phase shifts can then be obtained through L\, 

" kF{(ka) - LjF^ka) 



5i = tan 



-i 



LjGi(ka) - kG'iika) 



(29) 



where F and G are Coulomb WFs. 

The difference between our derivation and that of Canto and Brink is that we do not 
assume that the trial WF basis, ^ of eq.(|I|), can be factorized into the internal, relative 
and center of mass parts in the interaction region. In fact, it is impossible to do so for our 
delocalized WF. 

The main task remaining is to calculate the GCM kernel K GCM (s, if). The color part is 
standard. The spin-isospin part is also unaltered for the N-N channel, but different from 
that of others E3 for the N-A and N-E channels because we use a spin-flavor asymmetric 
hyperon WF and the antisymmetrization is restricted within the five u, d quarks. The 
orbital matrix elements are more involved due to derealization; in particular, the center of 
mass motion must be properly eliminated. We use a momentum projection method to 
project out the P c = part, 
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As mentioned before, our GCM basis WF includes not only l 3 r 3 , but also Z 6 r°, l 5 r x , l 4 r 2 , 
l 2 r A , / 1 r 5 , Z°r 6 configurations. Here l ni r n2 means, 



ni ni+112 

n - si) n (f) R {fj - s 2 ) 

i=l j'=ni+l 



which can be factorized to the form 



$Bi(6)$B 2 (6)exp<|-— 



n 



(f- s) 2 + n(r c - s c f 



(31) 



(32) 



where 



r = R 2 -R\, s = s 2 -si, 
_ m^i + n 2 i? 2 E"=i n 



7? 



ni 



i?9 



n 

j=m+i r i 
n 2 



nisi + n 2 s 2 
n 



(33) 



(34) 



and £i(£ 2 ) are t ne internal coordinates of the ni(n 2 ) quark cluster. 

Different configurations have different s c , but the same r c = (EiLiTi) Let us introduce 
the parameter center £ = (si + s 2 )/2. Then s c = i — (3 — i)s/6, where i — 0, 1, . . . , 6 
corresponds to the (ni, n 2 ) = (6, 0), (5, 1), (4, 2), (3, 3), (2, 4), (1, 5), (0, 6) particle partitions. 
(This result can obviously be extended to any N-particle system.) By means of t and s, the 
exponential part of eq.fl32|) can be written as 



exp 



f 

262 



nvriz,-, ^ 9 -., (3 — i) 2 9 n(3 — i),_ -, 
'f- s) 2 + n(f c - t) 2 + n- -s 2 + — -(f c - i) ■ 



n 



36 



(35) 



Due to the appearance of f c only in the combination (r c — i) in eq. (|35|) , the momentum 
projection, eq. (p0[) , can be written as 



df c e t * vrc tf(r 1 ,...,r 6 ) 



*(n,...,r 6 ; 



P c =0 



d^(ri,...,r 6 ) 

(36) 



The spurious center of mass motion part of the GCM kernel can thus be eliminated by a 
double momentum projection 



- J dtd?(${s)\H - E\V{s)). 



(37) 



When we used this momentum projection method to calculate the matrix elements of the 
kinetic energy and the Galilean noninvariant Darwin term of the one gluon exchange 
Fermi-Breit interaction, we obtained the same analytic formulas as those of Fujiwara B. 
Our variational method has been checked with Fujiwara's numerical results ® as well. 
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IV. RESULTS 



Initially, we carried out a unified, parameter free (i.e., all parameters are determined by 
single hadron properties and the color screening constant is taken from lattice QCD,) 
model calculation of the N-N, N-A and N-E interactions. This was done to check whether 
the quarks delocalize reasonably in the different flavor channels to give rise to qualitatively 
correct N-N, N-A and N-S effective interactions. These effective interactions are shown in 
figs. 1-3 (NN ST = 01,10; N A ST = Of, if, NE ST = Of , if, 0§, l|) and indeed 
qualitatively reproduce the phenomenological results. That is, the intermediate range 
attractions, usually assumed to be due to meson exchange, are reproduced by the quark 
derealization in the QDCSM. In the N-N channels, this model even gives rise to 
semi-quantitatively correct effective interactions; figs. 4-5 show the 1 So, 3 Si and 1 D2 N-N 
phase shift fits, plf 

Figs. 2-3 shows the results for N-T, and N-A channels. For N-A, the spin triplet state is 
somewhat more attractive than the spin singlet; channel coupling adds a bit more 
attraction to this state but leaves the spin singlet almost unchanged. We note that and 
a He both have spin zero ground states and spin one excited states. One might interpret 
this as evidence that the spin-singlet N-A interaction is more attractive than the spin 
triplet. However, the situation is not so transparent due to the complications of the 
interactions of the four bodies involved and, in addition, A-E° mixing effects. The spin one 
ground state of the deuteron is certainly an indication that the spin triplet N-N 
interaction is more attractive than the spin triplet, and we might reasonably expect (by 
flavor symmetry) that this should hold true for all octet-baryon combinations. However, 
due to the strong tensor interaction from pion exchange, the Nijmegen OBE model F [l 
nonetheless includes a more attractive spin singlet N-N central interaction. Furthermore, 
there is a paucity of direct data on scattering in the Y-N channels, and widely differing 
relative strengths for the central interaction are all consistent with both the available data 
and the nuclear states referred to above. The QDCSM, on the other hand, predicts that 
spin triplet N-N and N-A interactions are stronger than spin singlet ones. Clearly, which 
central interaction is stronger in each case merits additional study. 

For AT-E, we find the strongest attraction in the 13 = f 1 channel, while the 13 = §0, and 
fO channels both have a little weaker attraction (single channel case), and the 13 = |l 
channel is repulsive. Channel coupling has little effect on f 1, and pushes f from 
attractive to repulsive. These show that the Af-E potential is more strongly spin and/or 
isospin dependent than the N-A potentials, which have a little weaker dependence on spin. 
These results are in qualtitative agreement with the calculations of OBE models |l|J^| and 
hybrid quark model calculations ||, except that our attraction for the AT-E(fl) channel is 
too strong. 

To check if one can obtain a semi- quantitative fit of N-A and N-E scattering by fine tuning 
of the model parameters, two kinds of adjustment of the color screening parameter, k, have 
been made: The first keeps the color screening parameter k for the u, d quarks unchanged, 
i.e., k u = = l.llfm -1 , but allows k s for the s-quark to vary; the second one keeps 
Ku = Kd = k s = k but allows the value of k to vary for the N-A and N-E channels. 
Figs.6, 7 and 8 show phase shifts for the NA ST = Of , if and NE ST = Of , if and Of, if 
channels with k s = g k u . These are similar, but not identical to other hybrid quark model 
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results [ffiJ. Figs. 9, 10, 11, 12 and 13 show integral and differential scattering cross 
sections [ P2"| , |2"3"[ for NA,pE + and A qualitative fit is obtained even though our model 
phase shifts differ from others. This feature occurs in the meson exchange model as well; 
the Nijmegen models D and F [jT| also have quite different phase shifts from each other. 

m 



We have shown that 



Quantitatively, the N-A total cross section is fit quite well, 
channel coupling does not change the N-A interaction very much; therefore, this good fit 
will be maintained even after N-A and N-E channel couplings are taken into account. The 
p-E + cross section found in the model is larger than the experimental value p3f and this 
channel does not couple to any others. Hence, some fine-tuning may be needed. 
Conversely, the p-E~ cross section found in the model is smaller than the experimental 
value. J£J For that case, a channel coupling calculation is needed to determine if the fit 
can be improved by strong channel coupling effects. 



V. CONCLUSION 

The QDCSM has been used to calculate the effective N-N, N-A and N-S interactions. 
Linear confinement, with a color screening constant taking from lattice QCD, has been 
used in this calculation. All other model parameters are determined from the properties of 
baryons. This means that we have a parameter-free model calculation for nucleon-baryon 
(N-B) interations. We find that the quarks delocalize reasonably in the different flavor 
channels (10 altogether) to induce qualitatively correct, effective N-B interactions except 
that the N-A (1/2,1) and N-E (1/2,1) channels have attractions that are somewhat too 
strong. For the N-N channels, this model even gives semi-quantitatively correct phase 
shifts in the 1 5o, 3 Si, and 1 D2 partial waves. After fine tuning the color screening constant, 
we find it also gives semi-quantitatively correct scattering cross sections for the N-A and 
N-S channels. 

Several points need to be improved upon and to be checked further. 

• 1. So far, only single channel dynamical calculations have been done. The effect of 
dynamical channel coupling must be checked, especially for the N-E channel. 

• 2. Only central interactions have been included in this calculation. Non-central 
interactions need to be studied and higher partial wave scattering should be checked 
correspondingly. 

• 3. A better fit of the existing N-B scattering data better would support the QDCSM, 
but achieving such an improvement is not be the primary goal of this model 
calculation. In fact, nucleon-hyperon scattering data is so sparse that different meson 
exchange models fit the data perfectly well, and hybrid quark models can be made to 
fit the data as well, if one is willing to fine tune as has been done for meson exchange 
models. Since the fundamental strong interaction theory is certainly QCD, the N-B 
interactions should be an excellent area in which to study non-perturbative QCD. 
Meson-baryon and quark-gluon descriptions both are able to describe the N-B 
interactions if we are willing to include the whole hierachy of meson and baryon 
excited states and the totality of quark-gluon interaction diagrams @. However, we 
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may well ask which one is the most economical approach for including major 
non-perturbative QCD effects and for paving the way to new strong interaction 
physics such as exotic quark-gluon states, strangelets and so on. That the QDCSM 
obtains a qualitative, and in some cases, even a semi-quantitative, fit to N-B 
scattering and few nucleon systems @ in its very naive version, might be an 
indication that it includes a substantial component of the true physics. The QCD 
basis of its model Hamiltonian needs to be studied further. 

• 4. Nucleon spin structure studies show that the pure valence constituent quark 
model is only a first approximation. However, quark-antiquark excitation Fock 
components are certainly present in the ground state of the nucleon This fact 
should be taken into account in any quark model approach to hadron structure and 
hadron interaction studies; the QDCSM also needs to be elaborated to include it. 

This research is supported in part by the Department of Energy under contract 
W-7405-ENG-36 and in part by the NSF of China. 
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Effective potential in MeV vs. baryon separation in fm for N-N channels with 

confinement parameter value, k = 1.111/m -1 . 
Effective potential in MeV vs. baryon separation in fm for iV-E channels with 

confinement parameter value, k = 1.111/m -1 . 
Effective potential in MeV vs. baryon separation in fm for N-A channels with 

confinement parameter value, k = 1.111/m -1 . 
Phase shifts in degrees vs. center of mass energy in MeV for 1 Sq and 3 Si N-N 

channels compared with data §3. 
Phase shifts in degrees vs. center of mass energy in MeV for 1 D 2 N-N channel 

compared with data @. 
Scattering phase shifts in degrees vs. center of mass energy in MeV for N-A 

with confinement parameter value, k s = ~k. 
Scattering phase shifts in degrees vs. center of mass energy in MeV for iV-E, 

T = 1/2 channels with confinement parameter value, k s = ^k. 
Scattering phase shifts in degrees vs. center of mass energy in MeV for iV-E, 

T = 3/2 channels with confinement parameter value, k s = |ac. 
Cross section in mb vs. incident nucleon laboratory momentum in MeV/c for N-A 

scattering, compared with data il. 
Cross section in mb vs. incident nucleon laboratory momentum in MeV/c for 

scattering, compared with data @. 
Cross section in mb vs. incident nucleon laboratory momentum in MeV/c for 

scattering, compared with data c§. 
Differential cross section in mb vs. cosine of center of mass scattering angle 

at incident E - laboratory momentum of 160 MeV/c for p-E - scattering, 

compared with data @. 
Differential cross section in mb vs. cosine of center of mass scattering angle 

at incident E + laboratory momentum of 170 MeV/c for p-E + scattering, 

compared with data B. 
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Fig. 10 pL Scattering Cross-sections 
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Fig. 11 pZ Scattering Cross-sections 

(K = (4/9)K) 
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Fig. 12 pE Differential Cross-sections 
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Fig. 13 pZ Differential Cross-sections 

(K = (4/9)K) 
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Fig. 1 Effective potential 

(for NN with k=1 .1 1 1 frrf 1 ) 
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Fig. 2. Effective potential 

(for NEwith k=1.111 frrf 1 ) 
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Fig. 3. Effective potential 

(for NAwith k=1 .1 1 1 frrf 1 ) 
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Fig.4 N-N 1 S and 3 S 1 phase shifts 

(k=1 .1 1 1 fm" 1 ) 
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Fig.5 N-N 1 D 2 phase shifts 

(k=1 .1 1 1 frrf 1 ) 
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Fig. 6 NA Scattering Phase Shifts 
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Fig.7 NZ Scattering Phase Shifts 

(K = (4/9)K) 
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Fig.8 NZ Scattering Phase Shifts 
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Fig. 9 NA Scattering Cross-sections 

(K = (4/9)K) 
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